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Abstract. We report about high spatial resolution ob- 
servations made with the IRAM Plateau de Bure inter- 
ferometer of the 12 CO(J = 1 — > 0) emission from a 
1.1 x 1.1 kpc plane-of-sky field on a spiral arm of Messier 
81. With a beam of 5" « 90 pc, we identify 6 giant molec- 
ular cloud complexes with virial masses of w 10 6 Mq, in- 
cluding one associated with a giant HII region. The de- 
duced -/V(H2)/ico ratios are about 3 times larger on av- 
erage than those measured near the solar neighborhood, 
suggesting that the complexes are not self-gravitationally 
bound except, possibly, for the complex associated with 
the giant HII region; they could be the average of several 
clouds of mass a few 10 5 M Q and diameter <100pc. The 
linewidths are very narrow with respect to the measured 
sizes, so that the size-linewidth relation for M 81 clouds 
is very different from that in the Milky Way. The narrow 
linewidths imply smaller virial masses than for Galactic 
complexes of the same size, and this is consistent with the 
weaker CO emission from the GMCs in M 81. The low ve- 
locity dispersion suggests a lower mean volume density in 
the cloud and, possibly, a smaller scale height of the molec- 
ular gas than in Galactic clouds of the same size. Com- 
parison of the interferometer and single-dish line profiles 
indicates that, at most, 30% of the single-dish emission in 
this field is from a widespread distribution of small clouds, 
and thus the population of molecular clouds is rather dif- 
ferent from that in the Milky Way. 

The H2 surface density in M 81 is low: although the 
region studied here is one of the richer molecular regions 
in the disk, the molecular surface density is much smaller 
than the interarm regions of M 51 for example. The HI 
gas dominates and can explain most of the extinction seen 
at optical wavelengths in this field. In some other fields, 
the HI gas cannot explain the observed extinction, but 
previous lower resolution observations detected little or no 
CO there. The present high resolution observations imply 
that the molecular medium in M 81 differs from that in 
the Milky Way. 
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1. Introduction 

To understand the star formation activity at both large 
and small scales, one must know the distribution of the 
atomic and molecular gas on both scales. The importance 
of molecular clouds as star cradles and the easier study 
of rich molecular galaxies like M 51, made low CO emit- 
ting galaxies, like M 81, puzzling - as long as one relies 
on CO as a molecular gas mass indicator. However Sauty 



ct al. (1997) recently showed that a typical Sa to Sc spi- 
ral galaxy has three to four times less molecular gas than 
atomic hydrogen. This result is different from the study of 



Young & Knezek (1989) who find a higher H2/HI ratio in 



earlier type galaxies compared to late types. The discrep- 
ancy may be due to the selection's criteria as the galaxies 
in the sample of Young & Knezek are mostly FIR and op- 
tically bright objects and thus more probably CO-bright 
objects. 

In M 81 the 12 CO(J = 1 -> 0) emission has been 
mapped from the nucleus out to a radius of 7 '.5 and the in- 
ferred mass of molecular hydroge n is ~ 10 8 M© (Brouillet 
et al. 1991 ; Sage and Westpfahl 1991 ), assuming a stan- 
dard N(H2)/Ico rati o. The HI mass within the same ra- 
dius is 5 10 8 M (Hine 1984) which leads to a molecular to 



atomic gas ratio of 0.2, a bit less than the mean value (0.3) 
derived by Sauty et al. for Sab galaxies, and much less than 
the average value (2) derived by Young & Knezek. 

However even if M 81 does not appear anymore so 
atypical, the distribution of molecular gas with respect to 
the other tracers of star formation still raises questions. 

Reichen et al. ( |1994{ ) from the distribution of UV 
2000 A emission find that massive star formation appar- 
ently does not occur preferentially at the locations of the 
highest gas surface density. Furthermore in our CO sur- 
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veys (Brouillet et al. 1988, 1991) we found a peculiar lack 
of emission near the most luminous HII regions, which 
would confirm that they are efficient in dissociat ing th e 
molecular gas or in pushing it out (Kaufman et al. 1989a ). 
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But there are still a few HII regions where the visual ex- 
tinction is too high if only the atomic gas is taken into 
account and which must then have associated molecular 
gas. The dilution factor in the IRAM 30 m beam can ex- 
plain the observed low TV (H2) which would mean that the 
CO is highly concentrated on these HII regions. 

Other evidence that the CO measurements may under- 
estimate the molecular gas mass is provided by the dust 



lanes. Kaufman et al. (1989b) deduce, from measurements 
of the visual extinction and the HI column density, that 
the dust filaments in M 81 are composed mainly of molec- 
ular gas. Though CO emission appears to be associated 
with the dust regions on a large scale, we did not detect 
any strong emission with the 30 m towards two particu- 
lar dust lane po sition s of high optical depth (A„ > 16) 
(Brouillet et al. |l99l| ). The mass of molecular gas is at 
least a factor of 3 to 5 times lower than that inferred by 
Kaufman et al. ( |l989b[ ). This inconsistency suggests that: 
(i) the gas-to-dust ratio is different from the Galactic value 
- but it seems unlikely that it is as low as required to 
explain the discrepancy; or (ii) the -/V(H 2 )/ico ratio is 
different from the Galactic value; or (iii) the molecular 
clouds have a low excitation temperature T ex . 



Indeed Allen and Lequeux ( |1993[ ) have detected ex- 
tended molecular clouds associated with dark dust clouds 
in the inner disk of M 31. These molecular clouds have a 
very low excitation temperature and display much fainter 
CO emission than Galactic GMC's. In such a case the CO 
line intensity cannot so straightforwardly be used as an 
indicator of the molecular mass. 

High spatial resolution observations are needed to pre- 
cisely locate the molecular clouds in M 81 with respect 
to the dust filaments and HII regions. Furthermore as the 
CO emission is weak, we need the most sensitive interfer- 
ometer. 

From our previous 30 m observations we identified the 
best field in order to detect molecu lar co mplexes in M 
81. In the N5 field (Brouillet et al. |1991| ) CO emission 
follows the East- West dust lane and at the (-20", -10") 
position the CO profile is double peaked. The HII region 
138 (Kaufman et al. |1987| ) lies 30" farther to the West. 
This HII region is not representative of the HII regions 
with too much extinction; half of the observed HI in its 
direction is already sufficient to account for its extinc- 
tion. We observed a position on the dust lane, covering 
the (—20",— 10") 30m position and most of HII region 
138. 

In this article we present the 12 CO( J = 1 — » 0) obser- 
vations (Sect, ft) and then discuss the inferred N (H2)/ico 
ratio in M 81 (Sect. 0). In Sect. we comment on the as- 
sociation of the molecular complexes and the HII region. 

2. Observations 

CO ( J=l— 0) line observations at 115 GHz were carried 
out towards M 81 with the IRAM Plateau de Bure In- 



terferometer in August and September 1995. We used 
the compact CD configuration which led to a synthe- 
sized beam of 5"1 x 4 "8; this corresponds to 89 pc x 
84 pc at the adopted distance of 3.6 Mpc (Freedman et al. 
1994| ). The primary beam is 44" HPBW. The four anten- 
nas were equiped with SIS receivers and system tempera- 
tures ranged from 300 to 700 K. 1044+719 and 0954+658 
were used as phase and amplitude calibrators. 3C273, 
3C84, NRAO530 and 3C454.3 were used as bandpass and 
flux calibrators. The six units of the correlator were split 
as follows: two units covering 160 MHz each were used 
for the continuum, two 256 channels units and two 128 
channels units provided a velocity coverage of 155kms _1 
and 258kms~ 1 respectively, centered at V\ SI = 171 kms -1 
(V[sr = Vhelio + 6kms~ 1 ). The channel separation was 
0.156 MHz and 0.625 MHz respectively thus giving an ef- 
fective velocity resolution of 0.65 kms -1 and 2.6kms~ 1 
(since the channels are not independent). 

The phase center of the observations was a20oo — 
9 h 55 m 00?6, <5 2 ooo = 69°08'44"2. The coordinates were af- 
terwards transformed into the B1950 system and the phase 
center in this system is then aigso = 9 h 50 m 54!0, ^1950 = 
69°22'56"0. 

The total on-source integration time was about 34 
hours. No continuum emission was detected at a 5mJy 
level (3a). The line maps were cleaned using the CLARK 
algorithm and then corrected for primary beam at- 
tenuation. We achieved a sensitivity of ~0.5Kkms _1 
(Ict) at the center of the line integrated- intensity map 
(map integrated in velocity over a 52kms~ 1 range) and 
^2.5Kkms~ 1 (Ict) on the outside border. The flux has 
been converted from Jy/beam to K by multiplying by 
3.8K/Jy. 

We also made CO(J=l— 0) observations with the 
IRAM 30 m radiotelescope in November 1996. Due to the 
minimum spacing of the Plateau de Bure antennas (35 m), 
structures more extended than 15" cannot be detected. 
The 30 m telescope can thus perfectly provide the short 
spacings data to complete the interferometer set. We used 
two SIS receivers in parallel; the system temperature was 
about 500 K. Two 1 MHz filterbanks and an autocorrelator 
gave velocity resolutions of 2.6 and 1 .3 km s" 1 respectively. 
The half power beamwidth was 23". Unfortunately due to 
bad weather we could not complete our single-dish map 
in order to add it to the interferometer data. However, we 
obtained a map of size 20" x 20" with 8" spacing and a 
sensitivity of 3ct ^80mK (at the center) and 110 mK (on 
the edge) for 2.6kms _1 resolution. The pointing error is 
estimated to be 5" at most. 



3. Results 

Figure 1 displays the line integrated intensity map and 
Figure 2 the channel maps. The field, which is located 
in the northern spiral arm, displays a complex molecular 
structure and there is no clear velocity structure. 



N. Brouillet et al.: Identification of molecular complexes in M 81 



1 — i — i — i — i — I — i — i — I — r 
1 180.1 ' 



~i — i — i — i — I — i — i — i — i — I — i — r 
176.2 ' ' 




174.9 



69°23'20" 



69°23'00" 



69°22'40" 



Clean Beam 

J I I I I I I I I I L 



9 h 51 m 00 s 9 h 50 m 55 s 9 h 50 m 50 s 



Fig. 2. Channel maps of the CO (1-0) line emission before correction for primary-beam attenuation. The contour 
spacing is 40mJy/beam or 0.15K (the zero level is suppressed) and the clean beam is 5'/l x 4 '.'8. The cross refers to 
the phase center of the observations. The lsr center velocity of each map is written in the upper left corner and the 
velocity spacing is 1.3kms _1 . The rms noise before correction for primary-beam attenuation is 55 mK. 



From fig. @ we identified 6 molecular complexes, noted 
CI to C6 on fig. Ill and in Table pi They were defined as 
structures appearing on at least 2-3 adjacent channels. 
These complexes range in integrated CO flux from 1.3 
Jykms -1 for C6 to 5.0 Jykms -1 in the case of CI and 
C3, and in effective radius from 78 pc for C4 to 105 pc for 
CI. Complex CI is located near the phase center, but C3 



and C5 lie well beyond the half power point of the primary 
beam. 

In order to estimate the extended emission missed by 
the interferometer, we compared the 30 m spectra to the 
interferometric observations convolved to the same spa- 
tial resolution (see fig. g). The two data sets agree very 
well in velocity and line intensity except for the (16", 16") 
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Fig. 1. CO(l-O) line integrated-intensity map after correc- 
tion for primary-beam attenuation; the map is integrated 
in velocity from 145 to 197 km s -1 . The contour spacing is 
1 K kms -1 (the zero level is suppressed), corresponding to 
2 times the noise level at the center of the map, and the 
clean beam is 5'.'1 x 4 "8. The cross refers to the phase 
center of the observations at aigso = 9 h 50 m 54!0, chgso = 
69°22'56'.'0. The filled triangle marks the center of the HII 
region 138. 



position. As this position is at the border of the interfero- 
metric map and because of the noise, we cannot conclude 
if a GMC within the 23" beam of the 30 m was too far from 
the phase center to be detected by the interferometer or if 
extended emission is missing from the interferometer map 
at this peculiar position. Over the entire observed field 
the difference of flux is in the noise of the spectra and the 
missing flux can be estimated to 0-30% because of the 
calibration uncertainties. Thus there is no strong struc- 
ture more extended than 15" (300 pc) and no widespread 
distribution of small molecular clouds. 

We note that the lincwidths of the complexes are nar- 
row compared to Galactic complexes of the same size. The 
M 81 linewidths plotted on the velocity dispersion ver- 



sus size diagram of Solomon et al. ( |1987j ) are well below 
those for the Galaxy (fig. [|). For the structures found by 
the interferometer, R ~100pc, we should expect a width 
of 15-20 km s" 1 typically, while we find 7-14kms~ 1 . The 
lincwidths of the clouds in this paper (mostly 7kms _1 for 
a 5" beam) scale up to the lincwidths observed for M 81 
clouds by Brouillet et al. ( |!99lD with the 30 m 23" (400pc) 
and 12" beams (of the order of 5 to 20kms~ 1 ) via, approx- 
imately, the classical relation AV ex size 1 / 2 . This assumes 
the linewidths are due only to self-gravitation and turbu- 
lence of the molecular clouds (not to rotational gradients 
in the beam). In magnitud e the linewidths in both cases lie 
below the Solomon et al. ( 1987|) relation (or extrapolated 
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Fig. 3. Interferometric (top) and single-dish (bottom) 
maps convolved to the same spatial resolution (23"). The 
(0,0) position is the phase center of the interferometer 
data; the offsets are in arcsec. The velocity scale is 100 to 
250 kms -1 and the intensity T m b scale is -0.1 to 0.2 K. 



relation) for Milky Way clouds, e.g., with a 23" beam, we 
should have expected at least 30kms _1 for clouds similar 
to Galactic clouds, if the clouds filled the beam. But it 



is uncertain to extrapolate the Solomon et al. (1987) rela- 



tion to structures of 400 pc in size, the observed linewidths 
must include the velocity gradient across the beam, and 
if clouds do not fill all the beam, they will not show the 
complete V-gradient. 
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The M 81 linewidths are also narrow compared to those 
of the two molecular complexes identified by Viallefond et 
al. ( |1992D in M 33 and to the molecular complex mapped 
by Casoli et al. (1987) in M 31. On a larger scale, Table 
2 of Garcfa-Burillo et al. ( |1993| ) lists the CO linewidth- 
size relation for external galaxies and, except for NGC 
628, M 81 appears to have narrower lines. Although a 
weaker turbulence, perhaps resulting from a less active 
star formation throughout the M 81 disk, could explain 
the narrow lines, the star formation activity in M 81 (and 
NGC 628) is generally stronger than in M 31. Another 
possible explanation is the CO abundance: if the clouds 
are less optically thick, the lines can be narrower. 
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Fig. 4. Molecular cloud velocity dispersion a(v) as a func- 
tion of size S = 0.52R C as defined in Solomon et al. 
( |1987| ), see their figure 1. The dots represent the Galac- 
tic clouds and the corresponding fitted line is cr(v) = 
5°' 5 kms _1 pc _1 . The open squares represent the M 81 
complexes. For the M 81 complexes the uncertainty in size 
S is ±9pc, and the uncertainty in the velocity dispersion 
is ilkms" 1 . 



4. The N (H 2 )/J co ratio 

The value of the X = N (H 2 )//co ratio is still a subject 
of debate, especially because the "standard" value (X = 
2.310 20 H 2 cm- 2 /(Kkms- 1 ), Strong et al. |j~98g| ) adopted 
for the Milky Way does not seem to be va lid ev erywhere 
in the Galaxy. For example Digel et al. ( |l995| ) recently 
derived a lower value X = 1.06 10 20 H 2 cm _2 /(Kkms _1 ) 
towards Orion. The determination of X in other galaxies 
is thus of great interest. 

The resolution of molecular observations is increasing 
due to millimeter interferometers which allow us to resolve 
molecular complexes in nearby galaxies. The X values are 



greater towards very late- type or irregular galaxies. They 
also find that X depends on ico as X is larger for small 
Ico- Metallicity, in addition, is a fundamental parameter 
determining the X ratio, and is decreasing exponentially 
with radius in the Galaxy (Sakam oto |1996 ). From a study 
of 5 local group galaxies Wilson ( 1995| ) finds that X sig- 
nificantly increases as the metallicity of the host galaxy 
decreases. 

Our interferometric observations allow us to isolate 
molecular complexes in M 81. Though these structures 
are bigger than typical Galactic GMCs, we try to derive a 
value of X for these molecular complexes from two differ- 
ent methods: the virial theorem and the visual extinction 
Ay. The latter method requires an assumption about the 
relative distribution of stars and gas along line-of-sight. 

4--1. Virial mass 

In order to calculate the virial mass, we used M v - lr — 
95 D c AV 2 where M v - lr is in M Q , D c is the effective diame- 
ter in pc and AV is in kms^ 1 (Solomon et al. 1987). This 



not constant and Nakai & Kuno (1995) find the follow- 



ing tendency: compared to our Galaxy, the X values are 
smaller or comparable in galaxies earlier than Scd, and 



formula is obtained assuming that the cloud is a sphere 
with a power-law (a=l) density distribution. If instead we 
represent the cloud as a sphere with uniform density, the 
virial mass is then: Af v ; r = 104 D AV 2 and the mass is 
just 10% higher. 

To calculate the masses, we derived from the spectra 
the rms noise a and the linewidths at half intensity AV. 
And from the line contour maps we measured the size of 
the complex at a 3ct level as well as the integrated in- 
tensity Ico averaged over the cloud's area; an effective 
diameter D c is then derived from the area of the complex 

D 2 

(area=7r^p). We make the hypothesis that the clouds are 
spherical so that there is no effect due to the inclination 
of the galaxy. The cloud diameters were deconvolved from 
the synthesized beam using a gaussian model. Anyway as 
they are about twice the FWHM of the synthesized beam, 
the deconvolution just leads to a 10% smaller virial mass 
and a 10% smaller value of X. 

Because the clouds in this paper have a rather large 
sky-plane radius (~ 100 pc), it is likely that they are some- 
what flattened. The thickness of the molecular disk is 120 
pc (FWHM) in the Milky Way (Combes |l99lD and prob- 
ably smaller in M 81 because of the lower velocity dis- 
persion. Nevertheless, the following example suggests that 
the virial mass of the cloud calculated by assuming these 
clouds are spherical is close to that obtained for a some- 
what flattened cloud because the inclination correction 
tends to compensate for the oblateness correction. The 
virial mass for an oblate spheroid of uniform density and 
eccentricity e is the mass for the uniform density sphere 
multiplied by the factor e/arcsin e (Binney & Tremaine 
1987). The inclination correction to reduce the effective 
diameter to face-on (inclination i = 0°) increases from a 
factor of 1 for a sphere to a factor of (cos z) -1 / 2 for a 
highly flattened spheroid. The inclination i is 59° in M 
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81. As an example, if we suppose the clouds have a mi- 
nor/major axis ratio in the range 0.2-0.6, with the minor 
axis perpendicular to the plane of M 81, then, to get the 
virial mass of a uniform, oblate, spheroidal cloud in M 81 
from the virial mass of a uniform sphere, one multiplies by 
a total correction factor (for oblateness and inclination) of 
0.97-1.01, respectively. 

If we assume that the virial mass represents the molec- 
ular hydrogen and helium mass, we can derive the ratio 
X = N (H 2 )//co and the results are listed in column 12 
of Table 1. The total uncertainty on these values is about 
50%. 

The mean value of X is 6.7 ± 2.3 10 20 
H2cm _2 /(Kkms _1 ) for the six clouds. The uncer- 
tainties are much bigger for complexes C3 and C5 which 
are located on the border of the primary beam, and 
for complex C6 which has fainter emission. Omitting 
these clouds with less well determined parameters, the 
mean value of X is 6.1 10 20 H 2 cm- 2 /(Kkms~ 1 ). The 
mean value is about 3 times the standard Galactic 
value and about twice the va lue e xpected from the 
metallicity relation of Wilson ( 1995 ). The metallicity 
of this field is close to the one measured i n the solar 
neighborhood. Indeed Garnett & Shields ( 1993| ) find 
12+log(0/H)=8.94 towards HII region 138 and a mean 
value of 8.79 at the galactocentric distance of this field in 
M 81, whereas Shaver et al. Ql983| ) get 8.70±0.04 for the 
solar neighborhood. For such metallicities, as for M 31 



and inner clouds in M 33, Wilson ( |1995| ) find X = 3.3 10 20 
Hacm-VtKkms- 1 ). 

Except for cloud C3, the value of X is unusually large. 
If the molecular hydrogen + helium mass represents the 
virial mass, such a large ratio means that, though CO 
emission in this M 81 field is weak, there is a lot of H2. 
As the metallicity is normal, this could be due to a recent 
massive HI to H 2 conversion due to collisions; clouds CI, 
C2 and C5 are along the dust lane a little downstream 
from the HI ridge and may form a cloud collision front 
(Elmegreen 1988). We may be seeing only the warm en- 
velopes of cold molecular clumps, so the integrated CO 
intensity may not represent the entire molecular mass of 
the cloud. Alternatively, if the true X value in M 81 is 
in fact comparable to the Galactic one, the virial mass 
is about 3 times bigger than the molecular mass and ei- 
ther the clouds are virialized because of a large mass of 
stars and/or cold atomic gas not taken into account in 
these calculations (see next paragraphs), or the clouds are 
not virialized. In particular, Taylor & Wilson (1997) re- 
cently made interfcrometric observations with OVRO of 
another field on the spiral arms of M 81 with a higher 
angular resolution (3"). They have identified three GMCs 
with diameters ^100 pc and very narrow linewidths of 4.7- 
7.7 kms -1 . Their field also has nearly solar metallicity, 
and their values for the integrated CO flux and M v [ r im- 
ply that X is 1.4 ± 0.7 10 20 H 2 cm- 2 /(Kkms- 1 ) in their 
clouds. The GMCs' peak brightness temperatures are 3K 



and with beam dilution it is consistent with the mean 1 K 
peak brightness temperature of the complexes of this pa- 
per. Thus it is possible that the molecular complexes iden- 
tified in this paper contain several (2 or 3) gravitationally 
bound clouds of mass a few 10 5 M© and size <100pc not 
resolved in our data. In the case of cloud C3, which is as- 
sociated to the giant HII region (see sect. @), we find the 
standard X value and, though the derived parameters are 
more uncertain, it could be a single gravitationally bound 
cloud. 

In the virial mass we have neglected the mass of the 
stars. Indeed if we assume that the mass of the optical 
disk is about 10 11 M Q within a radius of 20kpc, with 
an exponential distribution for the surface density /1 = 
[ioexp(— r/ro), ro=3kpc, the surface density 6.5 kpc far 
from the center is ^ = 200M©pc~ 2 . This is comparable 
to a GMCs' typical mean surface density (170 M© pc~ 2 , 
Solomon et al. 1987), but the stellar distribution has a 



scale height of about 10 times that of the GMC, so that 
the volume density is about 10 times smaller for the stars 
than for the GMCs and we can neglect the mass of the 
stars at such a distance from the center. However the 
molecular surface densities derived for the M 81 clouds 
are smaller (30-150 M© pc~ 2 ) with the derived X ratio 
and even smaller (7-30 M©pc -2 ) with the standard X 
value. Thus the mass of the stars may be less negligible 
compared to the molecular mass and may contribute to 
the binding of the clouds. 

In these calculations we have also neglected the HI 
gas. However there is a good agreement between the CO 
and HI distribution (see fig. g) and velocities (see Table 
1, where the HI velocities were measured by Adler on the 
HI image with 12" resolution from Adler and Westpfahl 
1995). It seems that there is a shift between the HI and 



CO maxima but the difference of resolution (9" for the 
HI map of Hinc (1984) and 5" for the CO map) does not 
allow us to conclude about whether the shift is significant. 
In section 4.2 we conclude that for cloud CI the HI is not 
associated with the molecular cloud as it is mainly in a 
cold dust lane at high z; this is also true for cloud C2 
since the overlaid HI is along the same HI feature. In the 
case of cloud C3, it is clear that one should not include 
the entire HI column density in the cloud mass, because 
the resulting A/(HI) alone times the 1.36 correction factor 
for helium would equal M v i r (see sect. g). For the other 
complexes, if we take into account the HI correction in 
the virial mass calculation and we assume that all the 
HI column density measured towards the CO complexes 



(Hine 1984) also contributes to the virial mass M v i r = 
M(H 2 ) + M(He) + Af(HI), then we deduce a mean X 
ratio of 4.3 ± 1.8 10 20 H 2 cm- 2 /(Kkms" 1 ), or omitting 
complexes C3, C5, and C6, which have more uncertain 
parameters, the mean X is 3.0 10 20 H 2 cm _2 /(Kkms _1 ). 
Nevertheless, the most plausible explanation for the 
high values of X is that the complexes are not virialized 
but consist of a few unresolved self-gravitating clouds. 
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Table 1. Identified molecular complexes noted CI to C6 (see fig. [j]). The central coordinates are listed in Columns 
2 and 3. Columns 4, 5, 7 give, respectively, the rms noise, the central velocity, and the linewidth at half intensity 
calculated from the spectra averaged over the area of the clouds and with a velocity resolution of 2.6 km s" 1 . Column 6 
lists the central HI velocity (Adler, private communication). Column 8 gives the effective radius in arcseconds deduced 
from the area of the complex (area=7ri? 2 ) and deconvolved from the beam (1" corresponds to 17.45 pc). Column 9 
gives the peak brightness temperature measured from the spectra with a 2.6 km s -1 velocity resolution (1 Jy/beam = 
3.8K for the 5'.'1 x 4'.'8 synthesized beam). Column 10 lists the CO integrated intensity averaged over the clouds, 
column 11 the virial mass, and column 12 the deduced X ratio. 
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Fig. 5. HI contours (from Hine's map) superimposed 
on the CO line integrated- intensity map (grey scale). 
The first HI contour corresponds to a column density 
of 1.9 10 21 atom cm~ 2 (15M Q pc -2 ) and the spacing is 
0.3 10 21 atom cm -2 . The resolution of the HI map is 9". 
See fig. 1 for the CO map parameters. 



4.2. A v /T\g as 

Kaufman et al. ( 1989b| ) measured the extinction of sev- 
eral dust filaments in M 81 and compared them with 
the HI surface density. One of their positions located 
at a 1950 = 9 h 50 m 53?8A 950 = 69°23'00'/3 corresponds 
to our complex CI. To analyze their data they used a 
dust-cloud model where the cloud is located at a height 
z from the midplane and has a uniform distribution of 
dust. For CI they found z=200-400pc and r v =4-8. The 



HI colum n den sity AT (HI) is calculated from the HI map 
of Hine ( 1984 ) averaged over the VLA 9" beam and 
iV(HI) = (2.25 ± 0.75)10 21 atom cur 2 for CI. If we know 
the A v /N sas ratio, we can deduce A gas and then derive 
the molecular amount NCH2) necessary to account for the 
extinction. 

The value of N gas /E(B — V) varies from one galaxy to 
another because of the metallicity (Nakai & Kuno 1995), 
but the value of A V /E(B — V) should not vary with galax- 
ies if the dust properties are the same on a large scale and 
in a moderately dense medium. So we can use the Galac- 
tic gas to extinction ratio (Bohlin et al. [1978 ) corrected 
for metallicity (as in Kaufman et al. |1989a ) . We thus have 
A v /A gas =0.58 10~ 21 mag atom -1 cm 2 at i?=6 kpc. 

With A gas = 1.9 10 21 r v atom cm~ 2 , then using the 
HI column density measured with 9" resolution and as- 
suming all the gas is in front, we have A(H2)=A gas - 
A(HI)=5.210 21 atomcm- 2 for r v =4 or 1.3 10 22 atomcm" 2 
for t v =8. From our data we measured a CO integrated 
intensity 7co = 4Kkms" 1 . For this extinction compari- 
son, we give the integrated CO intensity centered on the 
dust position and averaged over the 9" FWHM of the HI 
beam, whereas Table 1 lists the integrated intensity aver- 
aged over the cloud. If we assume that the interferometer 
missed at most 30% of the CO flux, we get X in the range 
4.5 10 20 to 1510 20 H 2 cm- 2 /(Kkms- 1 ). If we average the 
CO emission over just the width of the dust lane, then 
^CO = 7Kkms~ 1 , and assuming that the interferometer 
detected all the emission in this small region, we get X — 
3.5-8.8 10 20 H2cm~ 2 /(Kkms _1 ). These various estimates 
are consistent with the values of X derived above. 

We assumed here that the molecular cloud is associ- 
ated with the dust lane and thus at a height z=200-400pc, 
which is unlikely. The molecular scale height in our Galaxy 
is about 50 pc and the star formation is more violent than 
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in M 81, so that the stellar winds are more present to push 
clouds up. 

An extinction of t v —A implies a gas column density of 
7.6 10 21 atomcm~ 2 . If the HI were concentrated in the lin- 
ear dust lane of width 2.5", the HI column density would 
be 8.1 10 21 atom cm -2 , enough to explain alone the extinc- 
tion. But the HI column density is similar towards neigh- 
boring regions, so there is no evidence for the HI to be 
confined in the dust lane. For cold and dense HI clouds 
with iV(HI)=810 21 atomcm" 2 , a typical Av=10kms _1 
and temperature T—10K, the HI is highly optically thick 
and we don't measure the true HI column density. We can 
finally assume that there is both diffuse optically thin HI 
(at a typical temperature of 100 K) and the HI confined 
to the dust lane at a high z and optically thick (at a typ- 
ical temperature of 10 K). It can be the same for the CO 
emission as we detect the less optically thick emission and 
we don't see the gas providing most of the extinction. 

We can thus account for the observed extinction given 
the measured HI and CO emissions. However, we note 
that the dust position in cloud CI is not one of the po- 
sitions measured by Kaufman et al. ( |1989b[ ) with a large 
discrepancy between the extinction and iV(HI). 

5. Cloud C3 and HII Region 138 



Cloud C3 and the Taylor & Wilson ( |1997D cloud MC2 
represent the first identifications in M 81 of CO emission 
associated with an individual giant HII region (see fig. 1). 
Previous CO observations of M 81 had too low a spatial 
resolution to connect the weak CO emission with the star- 
forming regions, so it has long been a puzzle as to how M 
81 is forming massive OB associations in the presence of so 
little molecular gas. The answer for the interferometer field 
of fig. 1 is that most of the molecular mass, as indexed by 
the CO emission, is in the form of 10 5 — 10 6 M complexes, 
suitable as birthplaces for large OB associations. 

The giant HII region 138, with excitation parameter 
U = 2 50pccm~ 2 and diameters 17"xl3" (Kaufman et al. 



1987]) appears inside cloud C3. In projection, the center of 
the cloud complex is displaced 7" = 120 pc from the cen- 
ter of the HII region. The location of the HII region is 
suggestive of a blister; possibly, the molecular gas on the 
eastern side of the HII region was photodissociated by the 
OB stars. Also, Kaufman et al. (1987) conjecture that a 
nonthermal radio source adjacent to the southeastern side 
of the HII region may be an SNR, which would have af- 
fected the molecular cloud. In the Milky Way, individual 
HII regions with free-free radio luminosity similar to that 
of HII region 138 are associated with molecular clouds of 
mass 410 5 to 410 6 M, 



1986). If cloud C3, 



i0 (Myers et al 

associated with HII region 138, is a single virialized cloud, 

its mass (0.910 6 M Q ) lies within this range. HII region 



3 10 5 M Q , i.e., at the low end of the Milky Way range for 
clouds related to such luminous HII regions. Integrating 
the entire HI column density over the solid angle of the 
HII region gives an HI mass M(HI) of 1.1 10 6 M Q for HII 
region 138 and 2.3 10 6 M for HII region 172 (Kaufman 



et al. 1989a), not corrected for helium. In both cases, the 



172 (Kauf man e t al. 1987) near cloud MC2 in the Taylor 



& Wilson ( 1997 ) field has a free- free radio luminosity sim- 
ilar to HII region 138 , but MC2 has a virial mass of only 



molecular cloud is displaced towards the western side of 
the HII region, and the atomic mass in the region exceeds 
the molecular mass. 

A layer of gas with column density 1.2 10 21 atomcm~ 2 
in front of the HII region would suffice to explain its ob- 
served extinction, A v =0.7±0.3mag. Averaged over the 
HII region, the CO integrated intensity is 1.9Kkms _1 , 
which corresponds to a total column density through the 
disk of 5 10 20 H 2 cm" 2 if we use the value of X = 2.7 10 20 
H2cm _2 /(Kkms _1 ) found for cloud C3, and the HI col- 
umn density through the disk is 2.4 10 21 atom cm -2 . Thus 
there is more than enough gas to account for the observed 
extinction, and the HII region is probably either at the 
midplane of the gas layer or slightly to the near side. 



6. Conclusion 

We have mapped and identified in M 81 six molecular com- 
plexes with virial masses ~ 1-3 10 6 M and radii ~ 100 pc 
in an optimally-chosen field of deprojected size 1.1 x 
2.1 kpc (the inclination of M 81 is 59°). Comparing the 
virial masses with the CO fluxes of these clouds, we find 
a mean value of the N(R 2 )/Ico ratio 3 times stronger 
than the standard Galactic value and 2 times more than 
the value expected for such metallicities. The most plau- 
sible explanation is that the complexes are the average of 
several gravitationally bound clouds not resolved in our 
data. A possible exception is cloud C3 which is associated 
to the giant HII region and, though the derived parame- 
ters are more uncertain, it could be a single gravitation- 
ally bound cloud. Severe beam dilution in the single-dish 
surveys (e.g., with the NRAO 12m telescope) appears to 
account, in part, for the weakness of the overall CO emis- 
sion from M 81 in such surveys. 

However we observed here an optimal field as it over- 
lapped one of the stronger fields in the CO single-dish 
surveys. Its CO emission may not be typical of M 81 as a 
whole: the sum of the virial masses of the six complexes is 
8.710 6 M Q , which gives an upper limit to the H 2 mass of 
6.4 10 6 M Q if helium is omitted. (On the other hand, if the 
standard value of X applies, the H 2 mass is 2.310 6 M Q , 
omitting helium.) The upper limit is ^6% of the total 
molecular mass of M 81. The corresponding HI mass of 
this field is 1.6510 7 M Q , which leads to a molecular to 
atomic gas ratio of 0.14 if we use the standard value of 
X and neglect a diffuse component, or up to 0.6 if we use 
the sum of the virial masses and take the upper limit on 
the diffuse component from the comparison of the 30 m 
and interferometer data. Furthermore this field contains a 
giant radio HII region coinciding with one of the molecu- 
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lar complexes and the other clouds are on the UV-bright 
arm, so that there is adequate heating of the gas. But the 
observed narrow linewidths have still to be explained. 

The upper limit to the average molecular column den- 
sity in this optimum field on the arm (obtained from the 
sum of the Virial masses) is approximately 6M pc -2 . 
This is much smaller than even the interarm surface den- 
sity o f mole cular gas in the CO-rich galaxy M 51 (see Rand 
et al. |1992[ ). 

The 3 GMCs in M81 detected by Taylor and Wilson 
( |1997| ) in another field of the spiral arms have CO inte- 
grated fluxes of 1.6-3.2 Jykms^ 1 , diameters of ~ 100 pc, 
and linewidths of 5-8 km s _1 . For comparison, our 6 com- 
plexes have CO integrated flux of 1.3-5.0 Jy kms -1 , diam- 
eters of 160-200 pc, and linewidths of 7-14kms~ 1 . Thus 
our complexes have somewhat similar 12 CO (1-0) fluxes 
but significantly larger areas. In both cases, the linewidths 
are narrow as compared to Galactic clouds of the same di- 
ameter. This has some interesting implications. Because 
of the low velocity dispersion, it is possible that the scale 
height of the molecular gas is smaller than in the Milky 
Way. However, if the complexes are not significantly flat- 
tened, then the mean volume density in the M81 clouds 
must be less than in Galactic clouds. This would mean 
that the clouds themselves are different from those in the 



Milky Way. Taylor and Wilson ( |1997| ) speculate that their 
clouds might be complexes consisting of several subclouds, 
like the Orion complex. This is one way in which the mean 
volume density could be reduced. 

From our interferometric and single-dish data, it seems 
that this field in M 81 has no fewer molecular complexes 
than other comparable galaxies but lacks a significant 
widespread distribution of small molecular clouds (in the 
limits of the calibration uncertainties, i.e. less than about 
30% of the total flux) . This lack of small molecular clouds 
could be due, in this particular field, to a collision front 



(see sect. 4.1) where the small clouds have agglomerated 
into big ones. Another possibility is that, as the molecular 
surface density in M 81 is low, the small diffuse clouds 
have less self- shielding against CO photodissociation than 
the GMCs. 

In the regions of M 81 with little UV emission (e.g. 
the interarm regions and parts of the central 4kpc HI 
hole), there is little heating of the gas, so most of the 
CO gas there may be too cold to have been detected in 
the surveys thus far. Allen et al. ( 1997 ) reach the same 
conclusion from their picture of HI as a product of star 
formation process, and they speculate that the gas in the 
main part of the disk of M 81 is largely in the form of 
molecular hydrogen with a low excitation temperature. It 
is still uncertain whether the 12 CO (1-0) emission from M 
81 as a whole can be used directly as a reliable tracer of 
the molecular gas. 

We conclude from the observations presented here 
that the molecular medium in M 81 differs from that in 



the Milky Way. 
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